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(Dated: December 7, 2007) 

Abstract 

Using a sample of tagged Df decays collected near the D*^D^ peak production energy in 
e + e~ collisions with the CLEO-c detector, we study the leptonic decay Df — > t + u t via the decay 
channel r + — > e + i^ e P r . We measure B{Df — > r + ^ r ) = (6.17±0.71 ±0.34)%, where the first error is 
statistical and the second systematic. Combining this result with our measurements of — > 
and Df -> t + v t (via r+ -> vr + P r ), we determine f Da = (274 ± 10 ± 5) MeV. 
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In the Standard Model (SM), the decay rate of a pseudoscalar meson Pq^ to a lepton 
neutrino pair £ + i>e is given by 



T(P Qq - fu e ) = G2AV f fl rn Q - qm ] 1 - ) , (1) 



where Gp is the Fermi coupling constant, Vq q is the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix element, vtiq^ is the mass of the meson, and is the mass of the charged lepton. 
Because no strong interactions are present in the leptonic final state t + V(., such decays 
provide a clean way to probe the complex, strong interactions that bind the quark and 
anti-quark within the initial-state meson. In these decays, strong interaction effects can be 
parametrized by a single quantity, fp, the pseudoscalar meson decay constant. In the case 
of the meson, fp> s describes the amplitude for the c- and s-quarks within the to have 
zero separation, a condition necessary for them to annihilate into the virtual W + boson that 
produces the t+vg pair. 

The experimental determination of decay constants is one of the most important tests of 
calculations involving non-perturbative QCD. Such calculations have been performed using 
various models [H, 0,(3,0] or using lattice QCD [1,0] (LQCD). Trustworthy QCD calculations 
within the I?-meson sector would enable the extraction of \ V t d\ from measurements of B° — B° 
mixing, and \V u b\ from (the very difficult 0, measurements of B + — > t + u t . Precision 
measurements of the decay constants fp, and fp s from charm meson decays are an attractive 
way to validate the QCD calculations used in the £>-meson sector. 

Physics beyond the SM might also affect leptonic decays of charmed mesons. Depending 
on the non-SM features, the ratio of T(D + — > t + Vf) f£iPt ~^ could be affected [§], as 



could the ratio T(D+ — > T + u T )/T(Df — > /i+z^J [lOj, Any of the individual widths might 
be increased or decreased. In particular, a two-Higgs doublet model [s| predicts a reduction 
in T{D+ -> i + u t ). 

Among the leptonic decays in the charm-quark sector, — > t^vg decays are the most 
accessible as they are Cabibbo favored (\V CS \ ~ 1). Furthermore, the large mass of the r 
lepton removes the helicity suppression that is present in the decays to lighter leptons. The 
existence of multiple neutrinos in the final state, however, makes experimental measurement 
of this decay challenging. 

In this Letter, we report the most precise measurement of the absolute branching fraction 
of the leptonic decay — > r + z/ r , from which we extract the decay constant fn a usingEqJTl 
We use a data sample of e + e _ — > D^DJ events collected by the CLEO-c detecto r | l2l . Il3l . 



14],ll5( at the center-of-mass (CM) energy 4170 MeV, near D^DJ peak production [16j. The 
data sample consists of an integrated luminosity of 298 pb _1 provided by the Cornell Electron 
Storage Ring (CESR). We have previously reported [171] measurements of — > and 
Df — > r + v T (via r + — > 7t + u t ) with these data. 



From the interaction point (IP) out, the CLEO-c detector [12j, ]13|, ]l4j, [15(| consists of a 
six-layer vertex drift chamber, a 47-layer central drift chamber, a ring-imaging Cherenkov 
detector (RICH), and a Csl electromagnetic calorimeter, all operating in a 1.0 T magnetic 
field provided by a superconducting solenoidal magnet. The detector provides acceptance 
of 93% of the full 4n solid angle for both charged particles and photons. Charged kaons 
and pions are identified based on information from the RICH detector and the specific 
ionization (dE/dx) measured by the drift chamber. Electron identification is based on a 
likelihood variable that combines the information from RICH detector, dE/dx, and the ratio 
of electromagnetic shower energy to track momentum (E/p). Background processes and the 
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efficiency of signal-event selection are estimated with a GEANT-based 



(MC) simulation program. Physics events are generated by EvtGen [19 



18] Monte Carlo 
, and final-state 

radiation (FSR) is modeled by the PHOTOS [20j program. The modeling of initial-state 
radiation (ISR) is based on cross sections for D*^DJ production at lower energies obtained 
from the CLEO-c energy scan [l6| near the CM energy where we collect the sample. 

The presence of two Df mesons in a Df^Df event allows us to define a single-tag (ST) 
sample in which a Df is reconstructed in a hadronic decay mode and a further double-tagged 
(DT) sub-sample in which an additional e ± is required as a signature of leptonic decay, the 

being the daughter of of the r . The D~ reconstructed in the ST sample can either be 
primary or secondary from D*~ — > D~j (or D*~ — > ir°Df). (We also use charge-conjugate 
Df decays for the tag; in this Letter, mention of a particular charge also implies use of the 
opposite one.) The ST yield can be expressed as n$T = 2iV£>sTesT, where N is the produced 
number of D* Df pairs, B$t is the branching fraction of hadronic modes used in the ST 
sample, and esT is the ST efficiency. 

Our double-tag (DT) sample is formed from events with only a single charged track, iden- 
tified as a positron, in addition to an ST. The yield can be expressed as wdt = 2iV£>sT£>sG e DT, 
where Bsc is the signal decay (SG) branching fraction, eoT is the efficiency of finding the 
ST and the SG in the same event. From the ST and DT yield expressions we obtain 
Bsc = (n B T/n ST ) x (6st/ ^dt) = (^DT/e)/n S T, where e (= gdt/ ^st) is the effective signal 
efficiency. Since eoT ~ cst x csg (where esc is the SG efficiency), Bsc is nearly independent 
of the uncertainties in esT- 

To minimize systematic uncertainties, we tag using three two-body hadronic decay modes 
with only charged particles in the final state. The three ST modes are Df — > 4>tt~ , 
Df — > K~K*°, and Df — > K$K~ . The — > tc + tt~ decay is reconstructed by combining 
oppositely charged tracks that originate from a common vertex and that have an invariant 



mass within ±12 MeV of the nominal mass [211 ] . We require the resonance decay to satisfy 



the following mass windows around the nominal mass [2l|: <\> — > K + K~ (±10 MeV) and 
K*° — > K + 7r~ (±75 MeV). We require the momenta of charged particles to be 100 MeV or 
greater to suppress the slow pion background from D*D* decays (through D* — > irD). We 
identify an ST by using the invariant mass of the tag M(D S ) and recoil mass against the 
tag M TCCOi i(D s ). The recoil mass is defined as M recoil (D s ) = [(E ee - E D f) 2 - |p ee - PdJ 2 ] 1/2 ' 
where (E ee , p ee ) is the net four-momentum of the e + e~ beam, taking the finite beam crossing 
angle into account; (£t> s ,PdJ is the four-momentum of the tag, with Ejj a computed from 
PD a and the nominal mass [2l| of the D s meson. We require the recoil mass to be within 55 
MeV of the D* mass j2l|. This loose window allows both primary and secondary D s tags to 
be selected. 

To estimate the backgrounds in our ST and DT yields from the wrong tag combinations, 
we use the tag invariant mass sidebands. We define the signal region as —20 MeV < 
AM(D S ) < ±20 MeV, and the sideband regions as -55 MeV < AM(D S ) < -35 MeV 
or ±35 MeV < AM(D S ) < ±55 MeV, where AM(D S ) = M(D S ) - m Da is the difference 
between the tag mass and the nominal mass. We fit the ST AM(D S ) distributions to the 
sum of double- Gaussian signal plus second-degree polynomial background functions to get 
the sideband scaling factor, and use that scaling factor for DT events also. The invariant 
mass distributions of tag candidates for each tag mode are shown in Fig. [TJ 

The DT event should have an ST, a single positron (p e > 200 MeV) with no other charged 
particles, and the net charge (Q ne t) of the event is required to be zero. These DT events 
will contain the sought-after Df — > r + v T (r + — > e + u e u T ) events, but also some backgrounds. 
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FIG. 1: The mass difference AM (D s ) = M(D S ) — mp s distributions in each tag mode. We fit 
the AM(D S ) distribution (open circle) to the sum (solid curve) of signal (double Gaussian) plus 
background (second degree polynomial, dashed curve) functions. 



The most effective discrimination variable that can separate signal from background events 
is the extra energy (E extr3i ) in the event, i.e., the total energy of the rest of the event. This 
quantity is computed using the neutral shower energy in the calorimeter, counting all neutral 
clusters consistent with being photons above 30 MeV; these showers must not be associated 
with any of the ST decay tracks or the signal positron. We obtain -E ex t m in the signal and 
sideband regions of AM(D S ). The sideband-subtracted E ext rs. distribution is used to obtain 
the DT yield. 

The -Eextra distribution obtained from data is compared to the MC expectation in Fig. [2J 
We have used the invariant mass sidebands, defined above, to subtract the combinatorial 
background. We expect that there will be a large peak between 100 MeV and 200 MeV from 
D* — > '-fDg decays (and from D* — > n°D s , 5.8% branching fraction |2lJ). Also, there will be 
some events at lower energy when the photon from D* decay escapes detection. 

After the AM(D S ) sideband subtraction, two significant components of background re- 
main. One is from Df — > K\e r v e . If the K\ deposits little or no energy in the calorimeter, 
this decay mode has an -E ex tra distribution very similar to the signal, peaking well below 
400 MeV. The second source, other semielectronic decays, rises smoothly with increasing 
-E'extra) U P to 1 GeV. Estimates of these backgrounds are also shown in Fig. [2J The optimal 
signal region in -E ex tra for DT yield extraction is predicted from an MC simulation study. 
Choosing -E ex tra less than 400 MeV 1 maximizes the signal significance. The number of 
non-peaking background events in the E extTa signal region is estimated from the number 
of events in the sideband region above 600 MeV scaled by the MC-determined ratio c& of 
the number of background events in the signal region, &W, to the number of events in the 
sideband region, \ft^ . The number of peaking background events due to the — ► K^e + u e 
decay is determined by using the expected number from MC simulation. The overall ex- 
pected number of background events in the -E ex tra signal region (6) is computed as follows: 
b = C(,6W (data) + b(K ( [e + h' e ) MC , where M h )(data) is the number of data events in the -E ex tra 



1 Note that with our chosen cut of £" C xtra < 400 MeV, we are including Dg — > t + v t ^ as signal. 
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FIG. 2: Distribution of E cxtTa after AM(D S ) sideband subtraction. Filled circles are from data 
and histograms are obtained from MC simulation. MC signal and the peaking background (Df — > 
K^e + v e ) components are normalized to our measured branching fractions. 



sideband region and b(K^e + v e )uc is the number of background events due to — > K'[e + h , e 
as estimated from our MC simulation. The branching fraction for Cabibbo-suppressed de- 
cay -D+ — > K^e + u e has not yet been measured. We determine this quantity by measuring 
B(D+ -> K° s e + u e ) = (0.14 ± 0.06 ± 0.01)% using a sample of 38548 D+ decays (more tag 
modes are used to increase statistics). 

The ST yield, AM(D S ) sideband scaling factor, DT yield with 400 MeV cut, and the 
number of estimated backgrounds events are summarized in Table HI We find nsT = 12947 ± 
150 and tzdt = 102 ± 12 integrated over all tag modes. 

The signal efficiency determined by MC simulation has been corrected for a few small 
differences between data and MC simulation. We weight the mode-by-mode signal efficiencies 
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TABLE I: Summary of ST yield (ust), ST mass sideband scaling factor (s), DT yield (?t,dt) with 
400 MeV cut, and the number of estimated backgrounds (6), where n s is the yield in the ST mass 
signal region, and n B is the yield in the ST mass sideband. 



Tag 


Mode 




n ST 


n ST 


s 


n S T 


n DT 


n DT 


b 


n D T 


DJ 


— > 4>7T~ 




5243 


391 


0.997 


4853.0 ±75.1 


49 





8.8 ±0.6 


40.2 ±7.0 


Ds 


-> K~ 




9020 


3661 


1.010 


5321.0 ± 112.8 


55 


3 


8.6 ±0.7 


43.4 ± 7.6 


Ds 






3499 


710 


1.022 


2773.1 ±65.0 


24 


2 


4.0 ±0.4 


18.0 ±5.1 



by the ST yields in each mode to determine e = (71.3 ± 0.4)% for the decay chain Df — > 
r + z/ T — > e + z/ e i/ r z/ r . Using £>(r + — ► e + u e v T ) = (17.84 ± 0.05)% [2l[ we obtain the leptonic 
decay branching fraction B(D+ — ► r + z/ T ) = (6.17 ± 0.71)%, where the error is statistical. 

The non-positron background in the signal positron sample is negligible (0.2%) due 
to the low probability (~ 0.1% per track) that hadrons (n + or K + ) are misidentified as 
e + . Uncertainty in these backgrounds produces a 0.2% uncertainty in the measurement of 
B(D+ — > r + u T ). The secondary positron backgrounds from charge symmetric processes, 
such as 7T° Dalitz decay (ir° — > e + e~7) and 7 conversion (7 — > e + e~), are assessed by mea- 
suring the wrong-sign signal electron in events with Q net = ±2. The uncertainty in the 
measurement from this source is estimated to be 0.9%. Uncertainties in efficiency due to 
the extra energy cut (1.8%), extra track veto (0.9%), and Q ne t = requirement (1.3%) are 
estimated using a sample in which both the and D~ in the event are tagged with any 
of the three hadronic ST modes. 

We considered five semileptonic decays, D+ — > 0e + z/ e , r]e + i , e , r]'e + i , e , K°e + i , e , and 
K*°e + i> e , as the major sources of background in the E extT3j signal region. The first two 
dominate the non-peaking background, and the fourth (with K° L ) dominates the peaking 
background. Uncertainty in the signal yield due to non-peaking background (0.5%) is as- 
sessed by varying the size of the dominant Cabibbo-favored semileptonic decays by the 
precision with which they are known [2l| . Imperfect knowledge of B(Df — > K°e + u e ) gives 
rise to a systematic uncertainty in our estimate of the amount of peaking background in the 
signal region. This uncertainty comprises two parts. We estimate the K\ showering sys- 
tematic uncertainty using -^(3770) events in which the D° has been fully reconstructed in a 
hadronic mode and the D° decays as D° —>■ K^Tt + Tt~ . When this uncertainty is combined in 
quadrature with the uncertainty in the determination of B(D+ — > Kge + u e ), the systematic 
uncertainty on B(D+ — ► t + u t ) is 4.5%. 

Other possible sources of systematic uncertainty include n ST (0.8%), tracking efficiency 
(0.3%), positron identification efficiency (1%), and FSR (1%). Combining all contributions 
in quadrature, the total systematic uncertainty in the branching fraction measurement is 
estimated to be 5.5%. 

In conclusion, using a sample of decays collected with the CLEO-c detector, we obtain 
a measurement of the absolute branching fraction, B(D+ — > t + u t ) = (6.17 ± 0.71 ± 0.34)%, 
where the first error is statistical and the second is systematic. This is the most precise 
measurement of this branching fraction and does not depend on measurements of other D s 
branching fractions for normalization. The decay constant fu g can be computed using Eq. [T] 



with known values [21( of Gp = 1.16637(1) x 10~ 5 GeV" 2 , \V CS \ = 0.9738 2 , m Ds = 1968.2(5) 



2 We assume \V CS \ = \ V u d\ and use the value given in Ref. 21 1. 
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MeV, m T = 1776.99io;26 MeV, and the lifetime of r Da = 500(7) x 10" 15 s (errors from 
these input parameters are negligible and ignored). We obtain fo a = (273 ± 16 ± 8) MeV. 
Combining with our previous decay constant determination [l7[ of fjj a = (274±13±7) MeV, 
we obtain fo s = (274 ± 10 ± 5) MeV. Our measured decay constant is consistent with the 
world average fn = (294 ±27) MeV j2l| and another recent measurement fo a = (283±17± 
7 ± 14) MeV [22]. These results are generally higher than recent LQCD calculations fo a = 
(249±3±16) MeV @ and f Da = (241±3) MeV @j. The predicted suppression @ that would 
be caused by a charged Higgs seems to be incompatible with experimental measurements 
combined with LQCD calculations. 

Combining with our previous measurement [l7[ of — > t + v t (r + — > 7r + iVK we obtain 



B(D+ — ► t + v t ) = (6.47 ± 0.61 ± 0.26)%. Using this with our measurement [17| of D 



u + v a i we obtain the branching fraction ratio ff^l^ r , Vt \ = 11.0±1.4±0.6. This is consistent 
with 9.72, the value predicted by the SM with lepton universality [lOj, 111], as given in Eq. [TJ 
We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the A. P. Sloan Foundation, 
the National Science Foundation, the U.S. Department of Energy, the Natural Sciences and 
Engineering Research Council of Canada, and the U.K. Science and Technology Facilities 
Council. 
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